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1 Abstract 

An improved unified formulation based on the effective field theory is introduced for a spin-1 Ising 
model with nearest neighbor interactions with arbitrary coordination number z. Present formu- 
lation is capable of calculating all the multi-spin correlations systematically in a representative 
manner, as well as its single site counterparts in the system and gives much better results for 
critical temperature, in comparison with the other works in the literature. The formulation can be 
easily used to various kinds of spin-1 Ising models, as long as the system contains only the nearest 
neighbor interactions as spin-spin interactions. Keywords: spin-1 Ising model ; Effective field 
theory ;Correlation functions; Bond diluted Ising system; Random field distributions, 
Crystal field dilution 

2 Introduction 

Spin-1 Blume-Capel (BC) model [T][2] is one of the most extensively studied models in statistical 
mechanics. Beside the various real magnetic systems, the model or some extensions can describe 
many physical systems such as ternary alloys and multicomponent fluids such as 3 He — 4 He 
mixtures [3]. Although the BC model is a very simple model, it exhibits many multicritical 
phenomena such as first (second) order order-disorder phase transitions, tricritcal points. On 
the other hand by inclusion of random distributions of crystal field and/or external longitudinal 
magnetic field or bond dilution of the BC model, makes the model more valuable. For example, 
BC model with random crystal field can describe 3 He — 4 He mixtures in a random media i.e. 
aerogel [3] . Taking into account these randomness effects will change phase diagrams of the model 
drastically, then model would exhibit richer multicritical phenomena. 

The spin-1 BC model was studied in different lattices by a variety of methods such as two 
spin cluster approximation [5], Bethe lattice approximation [BJ, series expansion method [7J[5], 
cluster variation method [SJ[TD] , Monte Carlo simulations [TTHTi] . renormalization group (T5HT7] 
and effective field theory (EFT) {T8H2(T| . In order to simplify the application and improve the 
results of EFT, Kaneyoshi introduced differential operator technique |211f25"] and Du introduced 
the expanded Bcthc-Pcicrls approximation (BPA) [26H28] . 

On the other hand there exist several works dealing with the BC model with random crystal 
field and bond dilution. Among these works, [2§H55] examined the random crystal field system 
with EFT. Some other applied techniques about this system are, mean field theory [MH4*0] . cluster 
variation method [JT] j Bethe lattice approximation [32] , finite cluster approximation [4"51T4"4] , Monte 
Carlo simulation (45], pair approximation [46] and renormalization group method |47j . Also the 
bond dilution problem in (transverse field) BC model was investigated with similar methods such 
as EFT fjBHBU] , finite cluster approximation [BTHBB] , Monte Carlo simulations [BBJ , series expansion 
|67) , two spin cluster approximation [68H70] and renormalization group method [71j . 

As seen in this very short literature about the spin-1 BC model and some of its extensions, 
the EFT method for solutions of this model has an important place in the literature. Thus, the 
aim of this work is to introduce an improved general effective field formulation for the spin-1 BC 
model and its variants, which is capable of calculating all multi site correlations of the system, in 
contrast to the other EFT formalisms in the related literature. 

A typical EFT method or finite cluster approximation start by choosing a finite representa- 
tive cluster and using exact spin identity (e.g. given in [72]) which gives the thermal average of 
magnetization. In the process of evaluating this identity in the cluster, multi site spin correlations 
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appear. Almost all of the works in the literature treat these multi site spin correlations by using a 
decoupling approximation (DA) which completely neglects them. Recently we have treated these 
multi site spin correlations in the spin-1 BC model [73H75] for the specific lattices and showed that 
taking into account these correlations considerably improves the results. Beside the improvement 
of the results, calculation of these correlations allows us to obtain some thermodynamic quantities 
more accurately (e.g internal energy, specific heat). 

The aim of this work is to obtain a general EFT formulation via treating the multi site spin 
correlations which is valid for spin-1 BC model and its variants for any arbitrary coordination 
number, as long as the system contains only nearest neighbor interactions as spin-spin interactions. 
For this purpose, the paper is organized as follows: In Sec. [3] we present this general formulation. 
As an application of our formulation we give the results of the pure, bond diluted and crystal field 
diluted spin-1 BC model on three dimensional lattices in Sec. 21 and finally Sec. [5] contains our 
conclusions. 



3 Model and Formulation 

The Hamiltonian of our system is given by 

u = - E J *J a i*i - E^ s * - E ff ^ c 1 ) 

<i.j> i i 

where Sj is the z component of the spin and it takes the values Sj = 0, ±1 for the spin-1 system, 
Jij = Jji > is the ferromagnetic exchange interaction between spins i and j, D± is the crystal 
field, Hi is the external longitudinal magnetic field at a lattice site i. The first summation in Eq. 
([T]) is over the nearest-neighbor pairs of spins and the other summations are over the all lattice 
sites. Jij, Di and Hi may be given with certain distributions or they can have the same values for 
all pairs/sites, i.e Jij = J,Di = D, Hi = H. 

We consider a lattice which has N identical spins arranged. We define a cluster on the lattice 
which consists of a central spin labeled so, and z perimeter spins being the nearest neighbors of 
the central spin. The nearest-neighbor spins are in an effective field produced by the outer spins, 
which can be determined by the condition that the thermal average of the central spin is equal to 
that of its nearest-neighbor spins. 

We start with spin identities for spin-1 Ising system [75] which are given by 

where Tri is the partial trace over the site i, (5 = 1/ (fe^T), ks is Boltzmann constant and T is the 
temperature, {/j} is any function as long as it is not a function of the site i, Hi is the part of the 
Hamiltonian which includes all contributions associated with the site i. While the inner average 
brackets stands for the thermal average, outer average brackets (with subscript r) stands for the 
random configurational average which is necessary for including random distribution effects. 

z 

Let us denote sf as Performing partial traces in Eqs. ^ and ([3]) with Ho = — sq^~~^ Jqs s s 

5=1 

—go-Do — sqHq for central spin (so) and with H$ = — s$ Jqssq — ssl — qsDs — sgHg for perimeter 
spin (ss) then with using differential operator technique [23j will yield 

«{/o}so» r = (^({fo}fl [1 + ss sinh ( J 0S V) + q s (cosh(J 05 V) - 1)]^ Fi(x)| a=0 (4) 



«{/o}<Zo)) r = / ({ fo}J[ [1 + s s sinh (J 05 V) + q 5 (cosh (J oa V) - 1)] ) ) F 2 (x)\ x=0 (5) 
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(({fs}ss)) r = (({ft} [1 + s sinh ( J os \7) + 9o (cosh (Jo* V) - l)]» r F 1 {x + 7 )| a=0 (6) 
(({/«}3j)) r = ({{fs} [1 + so sinh (J oi V) + q (cosh (J M V) - l)]» r P 2 (z + 7 )|x=o (7) 

where V = d/dx differential operator, 7 = (z — l)h and h is the effective field per spin. The effect 
of the differential operator on a function is given by 

exp (aV)Fj (x) = Fj (x + a) , j = l,2 (8) 

with any constant a. 

Eqs. (H])-© are our fundamental equations and the functions in the equations are given by via, 

r( n rr, 2sinh[/?(x + A)] 2 cosh[/3(x + Hi)} 

Gi(x, Di,Hi) = — — — , G 2 (x,Di,Hi) 



as 



2 cosh (a: + exp (-/3A)' 2 cosh [/? (a; + A)] + exp (-/3A 

(9) 

F j (x) = JdDidHiPniD^PHiH^Gjix.D^Hi), 3 = 1,2 (10) 

where Pq (A) and (A) are probability distributions of the crystal field and longitudinal mag- 
netic field respectively. 

We can obtain the central spin magnetization and the quadrupolar moment from Eqs. ((4]) and 
([5]) by letting {/o} = 1 and expanding right hand sides of them. Both expansions have multi site 
spin correlations like ({s\S2 ■ ■ ■ <7fc-i9fc)} r where k < z. The main problem here is to handle these 
correlations. These correlations are treated with DA in most of the works in literature due to the 
mathematical difficulties. In DA these correlations are treated according to 

«*i*2 ■ • ■ <fe-i3*» r = «*i» r «*2» P • ■ ■ «« fc _i» r {{q k )) r (11) 

i.e neglecting all the multi site spin correlations. But these correlations can be obtained by choosing 
suitable form of {fs} in Eqs. ([6]) or (J7|). This is the main point of this work. If we choose a proper 
form of {fs} in Eqs. ((6]) or ([7]), this time other multi site correlations appear on the right of sides of 
Eqs. ([6]) and (O, which include the central site. By the same way, these can be obtained from Eqs. 
(HI or ([5]) with a suitable {/o}- At the end of this process we can obtain a set of linear equations 
which includes these correlations as unknowns. For a given set of Hamiltonian parameters, we can 
obtain the effective field (which represents the effect of the outer (z — 1) spins) by usual condition 
that 

«*0» r = «**»r- (12) 

Once effective field h is determined for a given parameters, we can solve the system of linear 
equations and obtain the whole correlations of the system for given parameters. Since the effective 
field is very small in the vicinity of the phase transition point, the critical temperature can be 
obtained by letting 7 — > 0. 

In the formulation presented here, instead of deriving all correlations in the system, we construct 
representative correlations. Let us denote the correlation ({siS2 ■ ■ ■ s p q p +iq p -^2 ■ ■ ■ Qn)) as 

((4 P) 4"~ P) )) r = ({sis 2 ■ ■ ■ s p q p+1 q p+2 . . . q n )) r . (13) 



For instance, if we choose n = 2,p = 1 in Eq. (|13|) for the system with z = 3, the corre- 
lation (\4 1) 4 X /) = (( s i<?2)) r represents all the terms {{siq 2 )) r , ((si<7 3 )) r , ((s2<7i)) r , ((s 2 g3)} r , 
(( s 39i)) r j (( s 392)} r - This means that our solutions are under the assumption of that all these 
correlations are equal to each other. In general, we make an assumption that all n perimeter site 
correlations which includes p number of s$ and n — p number of qs are equal to each other and 
we represent them by (^Sg^qg* as given in Eq. (fT3|) . The same assumption holds for the 

correlations which include central site and n perimeter sites, i.e all correlations which include sq 
and p number of ss and n — p number of qs are equal to each other and all correlations which 
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include go and V number of s$ and n — p number of qs are equal to each other, and we represent 
them as ^s s^ ) ^"~ p) ^ and ^<7o4 P)< 25"~ P) )) , respectively. 

Deriving correlations simply based on choosing {/q} in Eqs. (j4j or (0, or {f$} in Eqs. ([6]) 
or (J7|) and chosen {/o} and {fs}, will change the operator in the average brackets on the right 
hand sides of them. For example, by the choice of {/o} = s k in Eq. (j4j) one applies Sk to the 
expression ((so)) r - Applying si to the expression ((so)) r will create the correlation ((soSi)) r - After 
then applying s 2 to the {{sqSi)) t will create {(SQS1S2)) r correlation, and so on. Thus, we get the 
correlation ((s si . . . s z )) with successive applications of to the correlation ((so s i • ■ • s k-i)) r 
[k = 1, 2, . . . , z) in z steps. 

Since our starting points of the derivation process are Eqs. (H])-© we can write any correlation 
in terms of the related operator which is applied to the related function, as in expressions Eqs. 
(SJ)- (0) . Since the operators on the right hand sides of Eqs. ((4]) and (|5|) are the same, we will use 
the same operator for getting the correlations e.g. ((so s i s 2)) r and ((qoSiS2)) r - The same reasoning 
also holds for the derivation of the correlations from Eqs. ([6]) and 0. Thus, instead of deriving all 
correlations separately in sequential orders, we will concentrate on the evolution of these operators 
in that sequential orders. Let us denote the operator on the right hand sides of Eqs. (Q| and (|5|) for 
{/o} = 1 as 0o,o and the operator on the right hand sides of Eqs. (JS]) and © with 5=1, {/i} = 1 
as $o.o- With these operators we can write our fundamental equations given in Eqs. ©-(CO) as 
follows: 

((so)) r = ((e m )) r F 1 (x)\ x=0 (14) 
(<si)} r = <($oo>) r iMx + 7)L=o (15) 



«<Zo» r = «e o)) r ^(x)U =0 (16) 
((?i)) r = (($oo)) r f2(x + 7)U=o (17) 

where the operators are given as 

z 

©oo = 11 I 1 + s <5 sinh ( J os V) + qs (cosh ( J 0(5 V) - 1)] (18) 

5=1 

$oo = [1 + so sinh ( J 05 V) + g (cosh (J 05 V) - 1)] . (19) 

Let us derive correlations in sequential orders given below. At each step of the derivation, the 
operators evolve in the sequences as given below: 

(( s l))r -> (( s l s 2)) r -> •■• H- ((s 1 S 2 ...S z )) r 
$00 -> $10 ->...->• $ z -l, 1 ' 

<(so)) r -> ((soSi)) r -> ((s sis 2 » r ... ((soSiS2.--s z )) r , 21 , 

000 -> ©10 -> ©20 ... -> zO 

((s 1 S 2 ...Sfc+l)) r -» ((siS 2 ...Sfegfc+l)) r ■•■ ~> ((siS2---Sfe-m+ig , fe-m+2---gfc+l)) r 

$fc0 -> $fcl -> ■■■ ~> $/cm 

(22) 

((s sis 2 . . . s k )) r -> ((s s 1 s 2 ...s k ^iq k )) r -> ... -> ((s sis 2 . . .Sfe- m g fe - m+ i . . .q fc )) r . 

©fcO ©fel ~> ■•• ~> ©fcm- 

(23) 

In Eqs. (|22[) and (|23j) we use the self spin identities of spin-1 Ising system during the derivation 
process, which are given as follows. 

Siqi = Si, SiSi = qi, qiqi = qi, i = 0,l,...z. (24) 

According to Eq. (|24[) . applying si to the correlation ((sis 2 . . . s;<ft+i, • • • qk)) r will produce the 
correlation ((sis 2 . . . si-iqiqi+x . . . qk)) r where I < fc, which can be obtained by applying the 
((<&k-i,k-i+i)) operator to the function F\ (x + 7) and taking the value of resultant expression at 
x = 0. 
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Sequential orders given in Eqs. ([2"U1) - (|2~3|) . expose the recurrence relations for the operators as: 



$fe,o = s fe+ i$fc_i,o, fc = l,2,...z-l 

Bfc, = Sfe©fe-i,o, k = l,2,...z , , 

$fc,m = Sfc_ m+ 2*fc,m-X) m = l,2, ...fc 

©fc,m = Sfc-TO+l©fc,m-l, 771=1,2, ...fc. 

In order to obtain these operators let us start with the operator 0oo which is just the operators 
on the right hand sides of Eqs. Q and ([5]) with {/o} = 1. By using the assumption about the 
correlations which is explained above (under the Eq. (|13jl). we can write this operator as 

z n 

9„o = £5>™p4 P) ^ p) (26) 

ri=0p=O 

where 

A np = (*)(")n sinh ( J o 5 V) [] (cosh(J 05 V) - 1). (27) 

However, obtaining G^o from Ofc—1,0 by using Eq. (|25[) requires the determination of any term 
Sksyqg" 1 ^ which will appear on the right hand side of &ko- This requirement is also valid for any 
process for obtaining the operators by using the recurrence relations given in Eq. (|25|) . At this 
stage we use self spin correlations for the spin-1 Ising system given in Eq. (f24|) . Hence there exist 
three possibilities for the term SkSg^qg 71 ^ according to Eq. (|24l) : 



(i+l) (m) 
3 5 IS 
S l )J m ) _ J j(+l)jro-l) 



k > I + m 

4 l+i; ?r _i; . «<fc<« + m (28) 

Now, one strategy is to obtain Q^o from Eq. (|25p starting by Eq. (|26p is making first few 
iterations and trying to capture the general expression for the QkQ- However we follow a slightly 
different strategy here. Each of the iterations in Eq. (|25|) takes the coefficient A np which is the 



multiplier of some s^gg" 1 ^ and place it to another term s g l ^ qg™ as a multiplier. We envision that 



the iterations in Eq. ()25[) will generate certain movement of the coefficients A np between the terms 
si qi as a multiplier of them. Then if we can track every coefficient during the each step of the 
iterations, we can determine the place of any coefficient for that iteration step. Thus we can write 
an expression for any operator in any iteration step. 

For instance let us take the coefficient Aqq and track this coefficient during the derivation in a 
sequential order given in Eq. (|2Tj) . 



A 00Sg q s -> A 00 s s q s ->...->■ A oSg q s ->...->■ Aoos^ q 5 (Z\)) 

Eq. (|29p implies that, the coefficient Aqq is placed as a multiplier of the term sf^qg on the fe.f/i 

step of the derivation, i.e the operator expression &ko has the term AooSgqg = A 00 siS2 ■ ..Sfe 
on the right hand side. If we know the locations of the all A np coefficients at the k t h step of the 
derivation we can get an expression for &ko- This is the general strategy which is used to obtain 
all the operators in Eq. (|2T|) . 

By inspection and taking into account Eq. (|28p. in order to get an expression for the operator 
0fco, we can give the locations of the coefficients of the k t h step of the derivation in an order given 
in Eq. (|2Tj) as follows. 

Aw^sfq^. (30) 



For n > 



Ano -> S (k) (o) , ~ ( 31 ) 

syqy , fc > n 
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For p 7^ and even 

s (p-k) q (n-p+k) ; fc < p/2 




, p/2<k<n (32) 

, k > n 
and for odd p 

r s (^),(«-f+ fc ) , fc <(p+i)/2 

4 fc ~ 1) ^ n "* +1) . b+l)/2<fc<n (33) 

1 , * > n 

Detailed derivation of Eqs. (f3"0)) - (j3"3")) is given in Section [SJ From Eqs. (f5U|) - (|33|) we can get an 

expression for the operator ft^o as 

ft - (n (1) 4- n (2 ^ 4- nPl.NJ , ( A + nWl „(*0 „("-*>) 



+ Y, D nktVr k+1) + E E A np st k) «t p+k) (34) 

n— fe n— k,p— 2k — 1 

where A; = 1, 2, . . . , z and the coefficients are given by 



k—l k—1 n" k—1 n' 

A 

np 



n— n— lp— 2 n—lp—1 



2k-l" 2k-2' 

D nk = E ^"fc = E (35) 

p=2 p=l 

The / sign on the upper limits of the sums in (|35[) indicate that, sum runs over the odd indexes 
while the II sign on the upper limits of the sums indicate that sum runs over the even indexes. 

In the sequential order given in Eq. (|23[) . we can obtain the operator 0fc rn after a few iterations 
by using Eq. (|28|) , and it is given by 

ft - ( n [1) -X- n<- 2 A e ( fe -™)„( m ) _l_ n (3) (fc-m+1) (m-1) 

Ufe, m — \ iJ h + u k j s s q s + u k s s q s 
+£ (A n0 + D%)st m) <lt k+m) 

7i— k 

Z Z 71 

, V" ( A _L n^A (k-m+l) (n-fc+m-1) V"^ V"^ . (p-fe-m) (n-p+fc+m) 

n—k n—k,p—2k 

where k = 1, 2, . . . , z; m = 1, 2, . . . k. 

From Eqs. (JSJ) and (|7|) we can see the form of the operator 4> o defined in Eqs. (fT5j) and (fTTf 

as: 

*oo = Bo + B lSo + B 2 q (37) 

where 

B = 1 

Si = sinh(JoaV) (38) 
B 2 = cosh(JoiV) - 1. 

If we use the first line of the Eq. ([23)1 by taking into account Eq. (|2"5)) , we can obtain the operator 
<f>ko after a few iterations as, 

^ko—tiQSg q s + nis s s q s + r> 2 q s s q s {6\)j 
where fc = 0, 1, . . . , z — 1. After then, by the same way, using the third line of the Eq. (|2"5j) we get 
•Pfcm — -D0S5 ?5 + Bis s s <7 5 + B 2 q s s q s (41) J 



G 



where k = 0, 1, . . . , z — 1 and to = 0,1, ... ,k. By using Eqs. (|24p and (fH))) we can obtain the 
operator So^fcm and (jo^/cm as 

rf* I D i D \ (k—m) (m) D (fc — m) (m) 

so^fcm = (B a + B 2 )s s y s >q\ '+Sigo4 X W1 i 

3- / r-> , 7-> \ (fc— (m) , i -> (k—m) (m) \ ) 

qo<5>km = (B + B 2 )q s\ q\ ' + B lSo s} 'q y s ' 

where fc = 0, 1, . . . ,z — 1, m = 0, 1, . . . , A; and this completes the process of derivation of operators. 
Note that the definition given in Eq. (|40[) covers the definition given in Eq. (|39[) with m = in it, 
but the same thing does not hold Eq. (|36|) for Eq. (|34|) i.e the operator O^o have to be calculated 
from Eq. (|34|) . 

Now, since the operators and functions are well defined, we can construct the correlations. 
But there are more than one way to constructing any correlation. Thus we have to determine an 
appropriate expressions from the equalities below, for deriving correlations. 

The number of z/2(z + 3) — 2 correlations of the type ((sg° "^q^^ (fc = 2, 3, . . . , z and 

to = 0,1,..., fc) can be obtained by starting from Eq. (fl"5|) or p7|) as 

{{s { t m \ [ ^)) r = «*fc-i,m» r *i (x + 7) U=o, m = 0, 1, . . . , k - 1 (42) 

(( s t m) 1s m) )) r = «*fc-i,m-i» r F 2 (x + 7 ) | a =o, m = 1, 2, . . . , k (43) 

respectively, where fc = 2, 3, . . . , z. In order to obtain the complete set of correlations of this type, 
we can not use only one of Eqs. (|42[) and (|4"3")l due to the restrictions on to values of them. 

Since the operator Qk,m defined in Eq. (|36|) can not have m = value, the correlations of the 



type \\ s ° s 5// are obtained separately from the correlations of the type ^sos^ "^g^" 1 ) 
(to 7^ 0). z number of correlations of the type //sqs^\\ (fc = 1,2, ... ,z) can be obtained by 



starting from Eqs. (jT4|) or (fT5|) and they are given respectively as 

(( S o4 k) )> r = «eM» r fiWUo (44) 

((so4 fc) )) r = «*o**-i,o» r i ; iO= + 7)L=o. (45) 

z/2(z + 1) number of correlations of the type (J^sqs^ m ^qg m ^\ (fc = l,2,...,z and m = 
1,2, ... ,k) can be obtained by starting from Eqs. (fT4|) or (flTf and they are given as 

((^t m) 1 { S m) )) r = «e fe , m » r ^i (*) U=o (46) 
((so S r m) <?i m) )) r = ((so$k-i,m-i)) r F 2 (x + 7) | a=0 (47) 

respectively. 

In order to get the complete set of correlations of the type // 8o8g°~ m ^ q^jS including m = 
we can also use 

((sosf- m) qt ) )) r = ((so$ k -i,m)) r Ft (x + 7 ) | x=0 , m = 0, 1, . . . , k - 1 (48) 

with Eq. 07}. For fc = 1, 2, . . . , z, either Eq. @SJ with ?n = and Eq. (g7J| with m = 1,2, ... A; 
or Eq. (|4T|) with to = fc and Eq. (|48p with to = 0, 1, . . . , fc — 1 will generate number of z/2(z + 3) 
correlations which are the complete set of correlations including sq. 

Number of z correlations of the type (('ft'' 5 ^)) = 1) 2, • • • , z) can be obtained by starting 
from Eqs. (|15j) or (|16|) and they are given respectively as 

((<Zo4 fe) )) r = ((90^-1,0))^! (x + 7) |x=o (49) 
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((qos^))^ {{e kt0 )) r F 2 (x)\ x=0 (50) 

Number of z/2(z + 1) correlations of the type ((qoSg° m ^qg m ^\ (k = 1,2, ... ,z and m = 
1, 2, . . . , fc) can be obtained by starting from Eqs. (fTBj) or (JTTJ) and they are given as 

((qost m) <lt ) )) r = m, m )) r F 2 (x) \ x=0 (51) 

((qosf- m) q ( r ) )) r = ((qo*k-i, m -i)) r F 3 (x + 7 ) \ x=0 (52) 

respectively. Again, the complete set of correlations of the type (^(qoS g k ^\\ including 
m = can be obtained from 

((qos { s k - m) q { s m) )) r = «gb**-i, m » P ^ (x + 7 ) | a=0 , m = 0, 1, . . . , k - 1 (53) 



with Eq. (|52|) . Due to the restrictions of m in Eq. ([52]) and Eq. ([53]). we can obtain the complete 
set of correlations which include go, from Eq. (|53l) with m = and (|5"2"j) with m = 1, 2, ... fc or (|52p 
with m = k and (|53[) with m = 0, 1, . . . , fc — 1. 

It is possible to obtain number of 3z/2(z + 3) — 2 correlations from Eqs. (|42 |) -([53 ]) and then to 
construct the system of linear equations (which has correlations and fundamental equalities given 
in ([T4"]) -(|17 |) . as equations) which has the number of 3z/2(z + 3) + 2 equations. However, in order to 
ensure the correctness of the correlation expressions given in Eqs. (f4"2"]) - ([53"]) we have to investigate 
the spin-1/2 limits of them. It is a well known fact that the spin-1 system behaves like spin- 1/2 
system for large positive values of the crystal field i.e. in the limit of D —> oo. This means that in 
this limit, si = becomes as an inaccessible state for the spins on the lattice sites. As a result of 
this, in this limit, all qi terms in expressions become 1. This criteria impose that Eq. (|46| can not 
be used for deriving that correlations. This fact is explained in Section [7] in detail. 

From Eqs. (flH) - (|17|) and the definition of the related operators given in Eqs. (J^SJ) and ([37]) . 
we can obtain the fundamental equations of the central and perimeter spin magnetizations and 
quadrupolar moments as 

««o»r = EEC« ((s^qt P) )) (54) 

n=0p=0 r 
z n 

(^o)),- = EE^ ) ((4 P) 4 n - p) )) r (55) 

« Sl )) r = + « So » r + K V ({ qo )) r (56) 

((qi)) r = K ( o 2) + K[ 2) {(s )) r + Rf } ({q )) r (57) 
where (i = 1,2 and j = 0, 1, 2) 



Cnp — (A np ) r Fj(x)\ x= o 



(58) 



Configurational averages in Eq. ([55)) can be evaluated by the given probability distribution of 
bonds Pj (Jij) and the functions Fi(x) are given by Eq. ([T0|) . coefficients A np and Bj are given by 
Eqs. dUJ) and ([3"g]t. respectively. 

We can obtain the correlations which has only perimeter sites, by using Eq. PU|) in Eqs. (|42p 
and (03]) with fc = 2, 3, . . . , z. From Eq. (g5J) with m = we have 



W)) r - < ] ((^t 1) qf ) )) r + 4 1] ((^r^)),,, + ((qos?-%P)) r (59) 
and from Eq. (|43p with m = 1, 2, . . . , fc we obtain 



4 fc - m) <T)) =4 2) ((st m) q { r 1} )) +K? ) ((sost m) qi m - 1) )) +4 2) (( qo st m) qt- 1) ] 

(60) ' 
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For correlations that include sq and perimeter sites, we can use operator expressions Eqs. (|34p 
and (dTJ, in Eqs. (05 and (g7]) respectively with fc = 1, 2, . . . , z. From Eq. ((33]) we get 



'MH 0) ))_ = (4 M) +4 M) ) ((4%f )) +4 1 ' 3 ' ((•W)) +E (<$ + ^ 4) ) ((4 fc) <T fc) 



I I T 



+ E^ 5) ((4 fc ^r fc+1) )) r + E E ^((4 p - fc) 9 r p+fc) )) r (6i) 

n= fe n— k.p— 2k — 1 

where = 1, 2, . . . , z and from Eq. (|47j) with m = 1, 2, . . . , k 



({sost m \t ) )) r = (4 2) + 4 2) )({s st m \^r 1) )) r + K{ 2) (( ? o4 fe - ro) 4 m - 1) )) r (62) 

Similarly, for the correlations that include qo and perimeter sites, we can use the operator 
expressions given in Eqs. (|3~4"|) and (|3SJ), in Eqs. (|5U|) and (f5T|) . respectively with k = 1, 2, . . . , z. 
From Eq. (JOT]) 



»W)) r = (4"' + 4 2 ' 2 ') ((« 0) )} r +4 2 ' 3 > ((4'-"<!? , » r +E + iS 4 ') ((« 

n—k n—k,p—2k—l 

and from Eq. (|51|) with m = 1, 2, . . . , fc 

rf" m) 4 m) )) = f4 2,1) +4 2 ' 2) ) ((4 fc - m) 4 m) )) +4 2 - 3) ((4 fe - m+1) 4 m - 1) 



(fe) (n-fc) 



E4 2 fe 5) ((4 fc - 1) 4- fc+1) )) +E E ((4 p - fc) 4 n - p+fe) )) ( 63 ) 



E(^o ) +4r ) )((4 fc - m) 4"- fe+m) 

n—k 

((4 fc - m+1) 9 r fc+m - 1, )) p + E E^ 



n — k n—k,p—2k 

where the coefficients 4 an< ^ 4fc are gi yen by (i = 1, 2) 



(64) 



4 J) = (^)/iW|x=o> J = 1,2,3 
4£ = (^2}/iWU=o, J =4,5 



(65) 



and this completes the derivation of correlations process. 

Fundamental equalities in Eqs. ([54 ]) - ([57]) and Eqs. ([Sl)]) -(f6"4 ]) constitute linear equation system 
which includes number of 3/2z(z + 3) + 2 equations of correlations for the system with coordination 
number z. By appropriately indexing the correlations, this system of equations can be written in a 
matrix form and can be solved numerically. If the system is defined by coordination number z with 
bond, crystal field and longitudinal magnetic field probability distributions (Pj (Jij), Pd (Di) and 
Ph (Hi) respectively), then the functions in coefficients can be calculated from Eqs. ((9]), (|T0|) . the 
coefficients in the operators from Eqs. (|27j) . (|35|) . ([38]) and finally the coefficients in the correlation 
equalities from Eqs. (|58j) . (|65j) with using this distribution functions. 

Instead of deriving correlations by using Eqs. (|32"]). (|4"3"|) , (gj), ([17]), ([50]) and (|5"T]): using Eqs. 
(|42|) . ([43]) . ([47]). gU), ([52]) and ([53]) for correlations and Eqs. ([54])- ([57]) for fundamental equalities 
will generate another set of linear equations which gives results that are identical to those obtained 
in [23], which is just the extended BPA of the spin- 1/2 system to the spin-1 system. 

For completeness of our work, we give the DA formulation corresponding to Eqs. ([54]) and (|55[) 
with the help of Eq. ([TT]) as 

z n 

« s o» r = EE c «p <( s °»" «9o»r p (66) 

n=0p=0 
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««>»r = EE C S (^0))r P (67) 

n=Op=0 

By solving the system of nonlinear equations given by Eqs. (|66[) . (|67[) with the coefBcients in (|58l) 
we get DA results. 

4 Results and Discussion 

Since the aim of our work is to develop a general EFT formulation which is capable of calculating 
multi site spin correlations, we give only few results in comparison with DA and EFT approxima- 
tions in the literature. Unless otherwise stated, we mean by DA results that the results obtained 
by the formulation given in Ref. [52] and the EFT results that obtained by the formulation given 
in Ref. [23]. The results of bond dilution problem on a honeycomb lattice (z = 3) and crystal field 
dilution problem can be found in our earlier works |74[|75) . respectively. 

4.1 Pure System With Zero Magnetic Field 

In this simplest spin-1 system, all probability distribution functions can be given by delta functions 
as Pj (J i:j ) = S (J i:j - J), P D (A) = 5{D l - D) and P H (Hi) = S (Hi). The functions of this 
system can be calculated from Eqs. ([9]) and (|T0| by using these probability distribution functions 
for crystal field and magnetic field and they are given as 

Fl(x) = 2sinh ^) F2ix)= 2cosh ^) (68) 

1{ ' 2 cosh (/3a;) +cxp(- / 8£>)' 1 ' 2 cosh (fix) + exp {-(3D) { ' 

Since Jy = J for all bonds then (|58)) can be written as 

C ® = (n)(l) t sinh ( JV )] P [ cosh ^V) - l] n ~ p F k (x)\ x=0 (69) 

where k = 1,2. We can write the hyperbolic operators in terms of the exponential operators in 
Eq. (|69[) and by making binomial expansion of them we can obtain: 

ca? -(;)(?) SSg ( " r ( ( -* l( - + * - 2 * - W^MI- 

(70) 

if we apply the differential operator to the function F k (x) according to Eq. (JHJ we can get 

as - ( ( SSg ( V ) ( ( ^y—mr + ^m 

(71) 

Other coefficients for this case can be obtained by the same way from Eq. (|58]l as 

K ( k} - F k ( 7 ) 

K[ k} = 1/2 [F k (J + 7 )-F k (- J + 7 )] (72) 
^ fc) = 1/2 [F fc (J + 7 ) + F fe (-J + 7 )] -f fc (7) 

where k = 1, 2. 

In order to construct the system of linear equations, all necessary coefficients can be obtained 
with the help of Eqs. (|68| . ([7T]) and (|72|). We can investigate all phase diagrams by solving the 
system by following the procedure given in Section [3] 

The phase diagrams in (D/J — ksT c / J) plane for honeycomb (z = 3) lattice can be seen in 
Fig. 14.11 in comparison with DA and EFT. As seen in Fig. 14.11 the three approximation gives 
qualitatively same phase diagrams but we should call attention to quantitative differences. While 
EFT gives slightly lower critical temperatures than DA, the present formulation gives lower critical 
temperatures than EFT for all crystal field values. The quantitative differences can be seen in Table 
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14. II for critical temperatures at D = and Table B~T1 for tricritical points (D t / J, ksTf/J), at which 
the second order phase transition line and first order transition line meets. The D t /J values are 
almost equal for three formulations, however ksTt/J values of DA and EFT are close to each other, 
while the present formulation gives lower values for the fcsTt/ J. This fact can be seen in Table 
14.11 The superiority of the present formulation can also be seen in the numerical values of critical 
temperatures in the limit D — > oo, which are just the critical temperatures of the corresponding 
spin- 1/2 system [75]. 

The phase diagrams in the (D/J — ksT c / J) plane for the various three dimensional lattices ( 
z = 6 simple cubic and z = 8 body centered lattices) can be seen in Fig. [2] in comparison with DA, 
as a limiting case of the bond diluted system, i.e. c = 1. Since the quantitative and qualitative 
relations of the phase diagrams of the EFT formulation and present formulation for these lattices 
are the same as z = 3 phase diagram, EFT diagrams have not been shown for this lattices. 




Figure 1: The phase diagram of the pure system for z = 3 in (D/ J — k B T c /J) plane. Solid and 
dotted curves correspond to the second and first order phase transitions, respectively and solid 
circles denote the tricritical points. 



Table 1: The critical temperatures of the 
pure system with zero magnetic field ob- 
tained by DA, EFT and the present work 
for the crystal field D/J = 0. 



Lattice 


DA 


EFT 


Present Work 


3 


1.519 


1.428 


1.302 


4 


2.187 


2.114 


1.952 


G 


3.519 


3.466 


3.265 


8 


4.849 


4.809 


4.587 


12 


7.512 


7.484 


7.241 



4.2 Bond Diluted System 

Let us treat the bond dilution problem for homogenous crystal fields (i.e Pd {Di) = 6 (Di ~ D) for 
all i) and zero magnetic field (i.e Hi = for all lattice sites). For this system we assume that the 
nearest-neighbor interactions are randomly distributed on the lattice sites according to probability 
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Tabic 2: The tricritical points (D t / J,k B T t / J) of the pure system 
with zero magnetic field obtained by DA, EFT and present work. 



Lattice 


DA 




EFT 




Present 


Work 




D t 


k B T t /J 


D t 


k B T t /J 


D t 


k B T t /J 


3 


1.427 


0.632 


1.424 


0.614 


1.452 


0.528 


4 


1.890 


0.949 


1.887 


0.946 


1.912 


0.846 


6 


2.816 


1.581 


2.814 


1.583 


2.841 


1.454 


8 


3.741 


2.220 


3.741 


2.217 


3.768 


2.056 


12 


5.594 


3.480 


5.593 


3.482 


5.621 


3.315 



distribution function 

Pj (J«) = (1 - c) 5 (J l3 ) + c<5 (J ti - J) (73) 

where c is the concentration of closed bonds and < c < 1. 

Since the crystal field and the magnetic field distributions are same as the pure system with 
zero magnetic field, the functions which are used for calculating coefficients for this system will be 
The configurational average in Eq. ((58)) . 



r {k) 



\ P [■ n 

p ) n / dj ° sPj ( j ° s ^ sinh ( j ° sv i ] / dj ° sPj ^ t cosh ( j ° 5v ) - ^ Ffe ( 



5=1 S=p+1" 

can be taken by using Eq. (f73| in Eq. (f74j) and it gives 



(74) 



i n / v p 



Ci nJ = \ t 1 )[csmh(JV)} p [ccosh(JV)-c} n - p F k (x)\ x=0 (75) 



where k = 1,2. 

Applying the same procedure between Eqs. (|69l) - ([7T|l to Eq. (|75]l gives 

=(:)(;) ( n ;")(;)(j) fjhi*^^^ 

(76) 



By the same way, the coefficients in Eq. (|58[) can be written for this system as 
K { k) = F fc ( 7 ) 

^ fc) = c/2[^(J + 7 )-F fc (-J + 7 )] (77) 
= c/2[ J F 1 fc (J + 7) + F fc (-J + 7)]-c J F 1 fe ( 7 ) 

where fc = 1, 2. 

As an example of the results obtained within the present formulation in comparison with the 
results of DA, let us investigate the phase diagram for simple cubic and body centered cubic lattices 
in (D/ J — ksT c / J) plane for the bond diluted system. The general effect of dilution of the bonds in 
spin-1 BC model is to decrease the critical temperatures. At a certain bond concentration value, 
namely bond percolation threshold (c*), ordered state is completely destroyed and the system 
cannot exhibit a ferromagnetic phase. Bond percolation threshold value depends on lattice type 
for a given Hamiltonian and for c < c* system shows no phase transition at all. 

It can be seen in Fig. [5] that, T c values of our formulation is lower than those obtained by 
DA for all values of c, as in the pure system with zero magnetic field. This means higher bond 
percolation threshold c* value for the present formulation with respect to the DA. This can be seen 
in Table 14.21 Since the T c values of the present formulation is closer to exact ones in comparison 
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with DA, it is expected that this is also true for the values of c*. In Tabic 22] we see that, the 
difference in the bond percolation threshold values obtained in the present work and DA and EFT 
gets smaller as the coordination number increases. 

For D — > oo it is obvious that the system behaves like spin- 1/2 system. This implies that, the 
phase diagrams in Fig. [2] are getting closed on the right side, after the bond percolation threshold 
value of the corresponding spin-1/2 system. In general, since the calculation of the correlations 
within the EFT gives closer T c or c* value in comparison with the exact results than DA or results 
of other EFT approximations, we can say that c value at which phase diagrams getting closed on 
the right side in Fig. [2](i.c corresponding S - 1/2 bond percolation threshold values) will be closer 
to exact ones also. This values can be found in our earlier work |76j . 

The order of the phase transition can be determined by examining the variation of the order 
parameter with temperature. As seen in Fig [31 we can conclude that all reentrant behaviors which 
emerge for the negative D/J values are of first order. In addition, magnetization behaviors for 
D/J= 2.0 for some selected c values can be seen in Fig. [3] 

The other superiority of the present formulation is computability of the some thermodynamic 
functions which are related to the correlations (e.g. internal energy and thus specific heat). The 
internal energy and specific heat of simple cubic lattice for D/J = 2.0 with some selected bond 
concentration values are plotted in Fig. 01 The behaviors of the internal energy and specific heat 
are more reasonable than DA results, since DA neglects all multi site spin correlations. Besides 
it gives zero internal energy after the critical temperature, which is physically impossible. One 
important point about the specific heat curves is, appearance of a hump as c decreases. While 
c decreases, peak of the specific heat gets smoother and gives place to a growing hump after the 
bond percolation threshold value c*. 




D/J 



CD 



z=8 
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* \ J -**^ e ^ 

O ^* 
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-4 -3 -2 -1 

D/J 
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Figure 2: The phase diagrams of the bond diluted system in the (D/J — kBT c /J) plane with 
selected values of the bond concentration c. Solid and dotted curves correspond to the second 
and first order phase transitions, respectively, while solid circles denote the tricritical points. For 
each c value, black curves are obtained by the DA and colored curves are obtained by the present 
formulation. 



4.3 Crystal Field Dilution Problem 

As a final investigation let us investigate the system with diluted crystal fields. The random crystal 
fields arc distributed to lattice sites according to given probability distribution function, 

P D (A) = (1 - p) 5 (A) + p8 (A - D) (78) 

The system is under zero magnetic field and there is no bond dilution, i.e A = and Jij = J for 
all This changes only the functions defined in Eq. (fT0|) . Using Pjj (Hi) = 6 (Hi) and Eq. (|78|) 
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Figure 3: The variation of the magnetization with the temperature for the bond diluted simple 
cubic lattice, with selected values of the bond concentration c. 
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Figure 4: The variation of the internal energy (U) and the specific heat (C) with the temperature 
for the bond diluted simple cubic lattice, with selected values of the bond concentration c. 



in Eq. (fTQ)) will give 

= 2(l-p)smh(f3x) 2psinh(/?x) 
U ' 2 cosh (fix) + 1 2 cosh (/3a;) + exp (-/?£>) 1 j 

= 2(1 -p) cosh (/fa) 2p cosh (fix) 

2[ ' 2 cosh (fix) + 1 2cosh(/3x) +exp(-/3L>) 1 J 

By using Eqs. ([75)1 and (1771) with c = 1 and Eqs. ([79)1 and ([50]) we can construct the system of 
linear equations. 

The distribution function given in Eq. (|78|) distributes crystal fields in p percentage of lattice 
sites D and remaining 1—p percentage of lattice sites zero. Because of this, p = phase diagrams 
in (D/J — ksTc/J) plane will be the line parallel to the D/J axis and the value of the fcsT c /J 
is just the critical temperature of the pure system (i.e. system with no external magnetic field, 
homogenous distributed D and J ) at D = 0. The phase diagrams in this plane will evolve towards 
the phase diagram of pure system in the same plane, when concentration value goes to p = 1. 
During this evolution, up to a certain p = p* value, phase diagrams will not intersect D/J axis 
on the left side of the (D/J — }cbT c /J) plane (i.e. region with negative D/J values). This means 
that, system is in the ordered phase at low temperatures for all D/J values for p < p* . 
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Tabic 3: The bond percolation threshold 
values (c* ) of the bond diluted system with 
zero magnetic field, obtained by DA, EFT 
and present formulation for different lat- 
tices. 



Lattice 


DA 


EFT 


Present Work 


3 


0.333 


0.347 


0.391 


4 


0.249 


0.256 


0.282 


G 


0.166 


0.168 


0.181 


8 


0.124 


0.126 


0.133 


12 


0.083 


0.084 


0.086 



The phase diagrams of the lattices with coordination numbers z = 6, 8 (simple cubic and body 
centered lattice respectively) can be seen in Fig. [5] in comparison with DA phase diagrams. Again 
we can see from Fig. [5] that the T c values of introduced formulation are lower than those obtained 
by DA. This makes the critical concentration value p* lower than that of DA. The p* values for 
different lattices can be seen in Table 14.31 



CD 




D/J 



CO 




D/J 



(a) 



(b) 



Figure 5: The phase diagrams of the crystal field diluted system in (D/J — ksT c /J) plane with 
selected values of the crystal field dilution parameter (p). Solid and dotted curves correspond to 
the second and first order phase transitions, respectively and solid circles denote the tricritical 
points. For each p value, black curves are obtained by the DA and colored curves are obtained by 
the present formulation. 



5 Conclusion 

In this work we have presented a general formulation for the spin-1 Ising system with BC model with 
nearest neighbor interactions. Within the introduced formulation we can calculate the multi site 
spin correlations in a representative maimer. This capability shows its superiority in the results 
of the critical temperatures and other critical values (e.g. bond percolation threshold values) 
and behavior of thermodynamic functions (e.g specific heat) as a function of the temperature. 
Our formulation can be applied to the systems with quenched disordered systems. These type of 
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Table 4: The critical crystal field concen- 
tration values of crystal field diluted sys- 
tem with zero magnetic held, obtained by 
DA, EFT and present formulation for dif- 
ferent lattices. 



Lattice 


DA 


EFT 


Present Work 


3 


0.484 


0.476 


0.379 


4 


0.604 


0.611 


0.544 


6 


0.729 


0.745 


0.718 


8 


0.794 


0.810 


0.802 


12 


0.861 


0.875 


0.892 



disorders in real materials can be extremely important since they can induce important macroscopic 
effects on the thermal and magnetic properties of the material. The systems corresponding to these 
systems cannot be solved exactly in most cases. Hence, some approximations must be done. But, 
including approximations in the calculations means that loosing accuracy of the results. This is 
a disadvantageous situation from the experimental point of view of the material science, although 
making approximations simplify the solutions. Thus, it is desired to improve approximations as 
much as possible to give more closer solutions to exact ones. 

As shown in Section 21 our formulation gives the best approximated solutions of BC model for 
the spin-1 system within the framework of EFT. Besides, it can give more reasonable behaviors 
of some thermodynamic functions (e.g. internal energy, specific heat) in comparison to the other 
EFT formulations. 

On the other hand we believe that, it is very important to develop a general formulation which is 
capable of applying to the related systems with arbitrary coordination numbers. In conclusion, we 
hope that the formulation and results obtained in this work may be beneficial form both theoretical 
and experimental point of view. 



6 Appendix A 

During the derivation process of correlations in a sequential order given in Eq. (|21[) we envision 
that each coefficient A np in the operators moves to another correlation as a multiplier of it at each 
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step of the derivation. Thus, derivation process generate a movement of the coefficients A np . We 
can track the 'path' all coefficients during this derivation process and determine the place of all 
coefficients at the kth step of this movement. 

According to Eq. (|28|) there are three possible 'ways' for a coefficient. Let us call movements 
in these three possible different ways as movements in regime 1, regime 2 and regime 3. In the 
regime 1, movement goes on according to the first line of the right hand side of Eq. (|28|) . Similarly 
in the regime 2 movement goes as the the second line of the right hand side of Eq. ([2"51) and the 
movement in regime 3 goes as third line. 

In Eq. (f29|). we give the path of the coefficient ^4oo from the beginning to the end of the 
movement. This movement always stays in the regime 1 since at each step, the step number k 
satisfies k > I + m, where indices l,m defines the term as Sg^qg™ 1 ^ which is the starting point of 
the kth step. 

Now we give the paths of the remaining coefficients according to this reasoning. At the begin- 
ning of the movement, the coefficients A rM (n ^ 0) are placed at the SgQg . For the first step 
(k = 1) since < k < n the movement starts in the regime 2. From Eq. (f2"5)l we can generate the 
path of these coefficients while the derivation process in a sequential order given in Eq. (f2"Tj) . The 
first n step of these coefficients in the regime 2 are given as follows 

W^^M-^-^W (81) 

We see from Eq. (f5Tj) that at the end of the first n step, the coefficient A n o arrives at s s n \ s QS> . 
After then, the movement goes on in the regime 1 since the step k — n starts to fulfill the condition 

k > n 

4^(0 s (« + x )(? (0) s (^(0) (82) 

Thus for the coefficient A n o, we can conclude that for the steps k < n it is at the place s s k \g n ^ 
(Eq. ([ST])), and for the steps k > n it is at the place s { s k) qf ] (Eq. ([52])) i.e we arrive Eq. pip . 

Remaining coefficients A np (n,p > 0) are placed at s s qg n ^ at the beginning. Since for the 
first step k = 1 < n is satisfied, then the coefficient A np starts its movement in the regime 3 and 
it stays in this regime until the step r which is p — r = r i.e. r = p/2 according to Eq. (|28|) . 

„(p)Jn-p) . (p-l) (n-p+l) Jp-r) (n-p+r) (fi o\ 

b S IS ~~* b S Is ■ ■ ■ ~ * b 5 'is \°°) 

If the index p, which is the number of ss in the term s^qg 71 p ^ is odd, then we have r = (p + 1) /2. 

Thus we can say that, the coefficient A np (n,p > 0) is at s s p ~ k ^ q s n ~ p+k \ for the end of the step 
k < p/2 (for even p) and for the end of the step k < (p+ l)/2 (for odd p) i.e. the first lines of the 
right hand sides of Eqs. (J22) and ([531) . 

At the end of the movement given in Eq. (|53")) the coefficient A np (n,p > 0) is at the place 

s< f #5™ P+V ^ ■ After than the movement goes in the second regime since the step number k 
becomes greater than p — r. The movement stays in the regime 2 for t steps as, 

„(p-r) Jn-p+r) (p-r+1) (n-p+r-1) Jp-r+t) (n-p+r-t) , RA \ 

b S y<5 b 8 ... —> A d - q s [0<i) 

Since in the regime 3 and regime 2 keeps the places of the coefficients s^qg" 1 ^ to be I + m 
constant, at a certain step, step number will catch I + m then the movement goes on a regime 1 
according to the first line of the right hand side of Eq. ([2"5]) . This step number satisfies r + t = n. 
Then while for even p this condition reads p/2 + t = n, for odd p reads (p + l)/2 + t = n. We 
find that for even p, t = n — p/2 and for odd p, t = n — (p + l)/2. Thus the movement given in 
Eq. (|84|) ends with Sg n \g 0S> for even p and s s n ^qg for odd p. We conclude from Eq. (|84|) that 
for even p the coefficient A np {n,p > 0) is at the place s s k \g n ^ and for the odd p at the place 

s s k qg n k+1 ^ for the end of the step r < k < r + 1. These results are the second lines of the right 
hand sides of Eqs. ([32]) and f33]) . 

Now the movement is in the regime 1 and according to Eq. ([2"5|) it goes as 
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From Eqs. and (|86|) wc can say that the coefficient A np (n,p > 0) is at the place Sg q$ for 
even p and at the place s$ ^Qg^ for odd p for the end of the step k > n. This corresponds to the 



third lines of the right hand sides of the Eqs. ([32]) and (|33|) and this completes our derivation of 
Eqs. 



7 Appendix B 

As explained in Section [3j Eqs. (|42j) - (|53j) have to give correct spin-1/2 equivalents in the limits 
D — > oo and — > 1 where i = 0, 1, . . . , z. 

The central spin and the perimeter spin averages for the nearest neighbor spin-1/2 Ising model 
with random bond and random magnetic field distribution can be given as |76j 

{{sv)) r = ((®o /2) )) r f(x)\x=o (87) 
({si)) r = ((^ 1/2) }) r /(^ + 7)U=o (88) 

where 

g(x, Hi) = tanh \p (x + Hi)], f(x) = J dHP H (H t ) g (x, H) (89) 
The operators in Eqs. (|87|) and (|88|) can be written as 

e< 1/2) = X>i 1/2) 4 n) (90) 
$ a/2) = B (i/2 )+jB (i/2) so (91) 

where the coefficients are given by 



A [ n /2) = ( M n sinh ( J 05V) J] cosh(J 0(5 V), n = 0,l,...z 

^ ' 5=1 S=n+1 

B { 1/2) = cosh(J 05 V) 
b{ 1/2} = sinh ( J 05 V) 



(92) 



With the operators 

$ (l/2) =jB (l/2) s ( fe)+jB (l/2) SoS W (93) 

S A m) -B^ So s^+B[^sf (94) 

ei 1/2) = 4 fe - 1) E4:/ + 2) 1 + 4 fc) E4I l /2) + t 4 1/2) 4- fc) m 

n=0 n=0 n=2fc+l 

we can obtain spin-1/2 multi spin correlations as 

«4 fe) » r = «C / 1 2) )) r /(- + 7)U=o (96) 

((-4 fe) », r -«er» r /WU=o (97) 

(( S °4 fe) )) r = ({»o*k-i)) r tt* + 7)l-o (98) 
Now, we can see from Eq. © 

lim Gi(x,#i) =g(x,Hi), lim G 2 (x,Hi) = 1 (99) 

D—>oc D—>oo 
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lim Fi(x) = f(x), lim F 2 (i) = 1 (100) 

D— >oo D— yoc 



and then 



where f(x),g(x) are the spin-1/2 functions which are defined in Eq. 

Let us label all operators and coefficients related to the spin-1 system with superscript (1). 
With a little combinatorics one can show that 

E^ffi - 4 1/2) ( 101 ) 

n—p 

is valid for the relation between the coefficients of spin-1/2 and spin-1 Ising system which is given 
in Eq. (I27D. 

Let us look at the qi — > 1, (i = 0, 1, . . . , z) limits of all necessary operators for deriving multi 
site correlations and fundamental equalities of spin-1 Ising system which are given in Eqs. (I20p , 
(|34f . (p6| . (|40j) and gj). We start with Eq. 061). 



hm eg = EE^4 P) = E4 P) E^ ( 102 ) 

n— Op— p—0 n—p 

According to Eq. (|TUT1) this will give Eq. (JOT]). Thus 

to8W = 8f (103) 
The spin-1/2 limit of the operator given in Eq. by using Eq. f[55j) can be obtained as 

: A;— 1 n" z 2k — l" 

e M = I E^io + EE4 5 + E E 4# ] 4 fc) 

=0 n=lp=2 n=k p=2 

k—1 n z 2k— 2 1 \ 2: n 

EE^y + E E 48 + E E A ^ [ r k) (m 

n=lp=l n=k p=l / n=k,p=2k— 1 

by rearranging the sums we get 

2k' z 2k— l' z z z 

hm eg) = 4 fc) EE48 + 4 fe - 1] E E48 + E E< )s ^ fe) ( 105 ) 

p—0n—p p—1 n—p p—2k+ln—p 

According to Eq. (|101j) . the right hand side of Eq. (|105[) is equal to the right hand side of Eq. 
(|95]l . then 

Kme«=ei 1/2) (ice) 

The spin-1/2 limit of the operator given in Eq. ([36)) is 

(z k—1 n" z 2k— l" 

E^io + EE48 + E E 48 ) 4 fe ~ m) 
n=0 n =lp=2 n=fc p=2 

'fc-1 n' z z 2k-2' \ z n 

EE48 + E<k- a + E E 48 4 fc ~ m+1) + E E 484 p - fc " m) dor) 

^n— lp— 1 n— A: n— A: p—1 / n—k,p=2k 

We can see from Eq. (|95|) that, 9fe„ m contains m , instead of sl fc m+1 ) i n Eq. (|107[) , then 
we can conclude that 

fimew^e^ (los) 
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The spin- 1/2 limit of Eq. (p7)) is given by 

*£>, - + b£ 15 ) 4 fc - m) + ^iV?- (109) 

With the help of Eq. (|38)1 , we can conclude from Eq. (|109[) with comparing that 

^^i = C^ (HO) 



Similarly we can obtain the results for the operators given in Eq. (|41j) as 

(1) _ , r *(V2) 



(1) _ ,(1/2) \ iLi ) 



Now we can look at the spin-1/2 limits of expressions given in Eqs. (I42l) - (I53[) . in order to decide 
whether they have correct limits or not. Let us start with Eq. (|42j) . Whereas the limit of the left 
hand side of Eq. (jUJ) is 

the limit of the right hand side of Eq. (|42|) can obtained as by using Eqs. (|100|) and (|110p 

^((€1,™)), ^ (* + 7) |«=o = (($i 1 _ / 1 2 L)) r / (» + 7) |x=o (H3) 
Then with the help of Eqs. (|112|l and (|113|l . we can obtain for (g2J 

:4 fe - W) )) =(M- / i 2 U) /(- + 7)U=o (H4) 



This is nothing but Eq. (|96|) for the correlation \\ s $ k "// s P m "l/^ svs t em - Thus we can 
conclude that expression given in Eq. (|42)) has correct spin-1/2 limit. 
In a similar manner, Eq. (|43|) has spin-1/2 limit as 

r m) » =((*£2)) i (us) 



From Eq. (|93|) we can write Eq. (|1 15|) as 



5 <5 

and from Eq. (1§2"1) . we have 



_// s (k-m)\\ B (^) 1 + // SQS f-m)\\ S (X/ 2)l (n6) 



(fc— m)\\ // (k—m) 



(117) 



i.e tautology. 

If we continue to work on the spin-1/2 limits of the expressions given in Eqs. (pH|) - (j53"|) . we can 
see that the expressions flU) , (gSJ) , (08]), (09J) and ((53j) gives the correct limits as in (|114[) for the 
expression (|42|) and the expressions (|47|) . (|50|). (|5T|) and ([52]) gives tautology as in (|117|) for the 
expression ([43]) . The only expression gives false spin-1/2 limit is (|46|). 
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